The quark-antiquark potential at finite temperature and the dimension two gluon 

condensate 
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A recently proposed phenomenological model, which includes non perturbative effects from di- 
mension two gluon condensates, is applied to analyze the available lattice data for the heavy quark 
free energy in the deconfined phase of quenched QCD. For large qq separations we recover previous 
results for the Polyakov loop, exhibiting unequivocal condensate contributions. For the qq potential 
at finite temperature and finite separation we find that a good overall description of the lattice data 
can be achieved once the condensate is properly accounted for. In addition, the model predicts a 
duality between the zero temperature potential as a function of the qq separation, on the one hand, 
and the quark selfenergy as a function of the temperature, on the other, which turns out to be 
satisfied to a high degree by the lattice data. 
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I. INTRODUCTION 

The understanding of the physics of the quark gluon 
plasma [l[ has triggered a lot of activity in the past years. 
(Seedi for accounts of the experimental situation.) As 
a part of this effort, there have been some recent advances 
in the study of the inter-quark forces at finite tempera- 
ture 0,[H|. In these analyses the change in free energy due 
to the presence of a static quark-antiquark pair separated 
by a distance r in a thermal bath has often been used. 
The color screening at high temperature will produce a 
dissolution of the quarkonium state, and this could be a 
signal of quark gluon plasma formation Q . 

The vacuum expectation value of the Polyakov loop is 
related to the propagator of a static quark and consti- 
tutes a natural order parameter for the deconfinement 
phase transition of quenched QCD. Due to ultraviolet 
divergences in the heavy quark selfenergy, the Polyakov 
loop is subject to multiplicative renormalization (see e.g. 
@> HI)- The perturbative computation of the renormal- 
ized Polyakov loop has been carried out to next to leading 
order [9( and only recently a practical systematic proce- 
dure has been devised to obtain it non perturbatively 
in lattice QCD, based on the computation of singlet and 
octet correlation functions (Wilson lines) • These cor- 
relation functions are related to the free energy of a heavy 
quark-antiquark pair. In the limit of large separation the 
two Polyakov loops decouple, and the result is essentially 
the square of one single Polyakov loop. In the opposite 
limit of separations much smaller than the Debye length 
the result is temperature independent, so that the qq free 
energy reproduces the known QCD static qq potential at 
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zero temperature, V q q(r), [111 ] . There is an additive am- 
biguity in the free energy that is fixed at short distances 
by comparison with this zero temperature quantity. This 
strongly suggests that there exists a deep connection be- 
tween the zero temperature qq potential and the vacuum 
expectation value of the Polyakov loop. 

The perturbative computation of the qq potential has 
been studied for a long time. At tree-level, V qq -(r) is a 
Coulomb potential arising from one-gluon exchange. Up 
to now there are estimates of the three loop perturba- 
tive contribution using several methods [l2|, [l3|, [l4| . The 
perturbative series turned out to be very poorly conver- 
gent at r > 0.1 fm and resummation prescriptions lead 
to large uncertainty due to scheme dependence. Several 
techniques have been proposed to improve the perturba- 
tive series: extraction of the C(Aqcd) infrared renor- 
malons [pH Urij . Operator Product Expansion (OPE) 
for r <C Aqcd -1 [1Z|, e ^ c - Non-perturbative methods, 
such as lattice simulations, and phenomenological poten- 
tial model analyses of experimental data for the heavy 
quarkonium spectra, indicate that the qq potential can 
be well approximated by the sum of a Coulomb term 
plus a linear term that becomes dominant at separations 
of the order ~ 1/A QCD PS El. This non-perturbative 
linear behavior at large distances is consistent with the 
quark confinement picture based on the QCD string. For 
a recent attempt within the holographic QCD approach 
based on the Ads/CFT correspondence see e.g. [2fJ and 
references therein. 

On the other hand, at temperatures of the order of 
the deconfinement phase transition the characteristic Eu- 
clidean scale corresponds to /i ~ 47rT c ~ 3 GeV, and 
one expects the OPE formalism to apply [21]. This ap- 
proach consists of expanding the Green's functions in 
inverse power series of p 2 , each term corresponding to 
a local operator, implying that condensates and power 
corrections could play a role. This has been verified 
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in a recent analysis of the renormalized Polyakov lat- 
tice data [1 01 ] . where the existence of non perturbative 
contributions driven by gluon condensates has been ex- 
posed (22| |. The first BRST invariant condensate is of 
dimension two {A^ a ), with the gluon field |23j. It 
is tempting to develop a model based on the computa- 
tion of correlations of Polyakov loops to connect these 
two regimes (zero temperature and large distances), so 
that the non-perturbative contributions in the Polyakov 
loop have the same origin as those in the zero tempera- 
ture qq potential, namely, the gluon condensates. This 
will be one of the issues of the manuscript. More gener- 
ally, we extend the one-gluon exchange model introduced 
in [22| for the Polyakov loop, to the description of the 
singlet heavy quark-antiquark free energy and compare 
with available lattice data. A preliminary study in this 
direction was presented in [24j ]. 

The paper is organized as follows. In section |TT] we 
review the observed non-perturbative behavior of the 
renormalized Polyakov loop just above the deconfinement 
phase transition, and the importance of thermal power 
corrections to describe it, thereby justifying the intro- 
duction of dimension two condensates in the theory. In 
section [TTT] we study the influence of these condensates 
on the free energy of the static qq pair. In section HVl we 
compare our model with available lattice data and estab- 
lish a remarkable duality between the zero temperature 
qq potential and the Polyakov loop. Finally, in section fVl 
we summarize our points and draw our main conclusions. 



II. THE POLYAKOV LOOP 

In this section we summarize results derived for the 
Polyakov loop and the role played by the dimension two 
condensate. Fuller details are given in [22j . The Polyakov 
loop or thermal Wilson line is the propagator of a static 
color source. Physically, the expectation value of the n- 
point Polyakov loop correlation functions are related to 
the change in free energy arising from the presence of 
static quark and antiquark sources in the heat bath. The 
Polyakov loop is defined by 



fl(x) = Texp ig 



l/T 



dx A (x,x ) 



(1) 



where T denotes time ordering, and xq is the Eu- 
clidean time coordinate. Aq = J2 a Ao, a T a and T a , a — 
1, . . . , AT? — 1, are the generators of SU(A r c ) in the fun- 
damental representation, with normalization tr(T a Tb) = 
(5 a h/2, and N c the number of colors. Let 



L(T) 



1 



-tin 



(2) 



be the normalized expectation value of the Polyakov loop 
(tr is the trace in the fundamental representation of the 
color group) . This quantity was computed perturbatively 



to next-to-leading order (NLO) by Gava and Jengo 0]. 
To leading order (LO), 



L(T) = 1 



1 

167T N c 



g 2 m D + 0(g 4 ) 



(3) 



where we have defined the dimensionless quantity rho = 
mo/T, and rau is the Debye mass. To one loop [25[ 

m D = ^=g(N c /3 + N f /6) 1 / 2 (4) 

Nf being the number of flavors. The temperature depen- 
dence in the perturbative L comes from the running of 
the coupling constant. 

In a static gauge (i.e., one where Aq is time indepen- 
dent), the Polyakov loop takes the simpler form 



fi(as) 



_ p igA ( x )/T 



(5) 



and the computation of L(T) can be undertaken using 
dimensional reduction ideas [S, S H3, S HI • Namely, 
after the non stationary Matsubara modes are integrated 
out, the static modes A (x) and A(x) are described by an 
effective three dimensional Lagrangian with computable 
parameters. 1 Expanding the exponential one finds 2 



^M4)> 



(6) 

The quartic term is 0(g 6 ) and hence, to order 0(g 5 ) one 
can use the Gaussian-like approximation 



L = exp 



9 2 (Ala) 

4N C T 2 



0(g 6 ). 



(7) 



This formula becomes exact in the large N c limit [30]. 
From here it is immediate to relate the Polyakov loop to 
the gluon propagator of the dimensionally reduced theory 



(Ao ta (x)A 0jb (y)) = S ab T 



d 3 k 



(2tt) 3 

The lowest order perturbative term 
£ P o(fc) ' 



ik < x -y^D w {k). (8) 



At 



(9) 



1 A subtle point is that to make contact with the Polyakov loop, 
the integration of non static Matsubara modes should be done 
in a static gauge, whereas most calculations of the dimensionally 
reduced Lagrangian are carried out in a covariant gauge. Fortu- 
nately the spurious gauge dependence so generated only affects 
higher perturbative orders 123] ■ 

2 Odd orders vanish assuming that the conjugation symmetry 
Afi(x) — * —A^{x) is not spontaneously broken. By the same 
token (Q) and (O t ) are real (and equal in the absence of bary- 
onic chemical potential, which we assume throughout). 
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FIG. 1: The logarithm of the renormalized Polyakov loop 
in gluodynamics. Lattice data from [l(J for A*',- = 4 (white 
squares) and N T = 8 (black squares). The straight line uses 
Eq. (|lip with the continuum extrapolated values of a and b 
in HU) (rather than the best fits to the N T = 4 or N T = 8 
data). LO and NLO perturbative calculations [j| are shown 
for comparison. 



(the upperscript P refers to perturbative contribution) 
yields 



(AL) P = -(N? 



i)^t 2 

; 4vr 



0(g 2 ) 



(10) 



where dimensional regularization rules have been applied 
to deal with the UV divergence at coincident points. 
When combined with the Gaussian-like formula (JT]) this 
term reproduces to LO the perturbative result in (|3j). 
The NLO result is also reproduced in this way (22J. 

Because in perturbation theory the dimensionless ratio 
T/Aqcd appears only through logarithms from radiative 
corrections, the finite temperature perturbative conden- 
sate (Aq a ) p is roughly proportional to T 2 and the per- 
turbative L(T) is rather flat as a function of T. This is 
definitively not what is seen in the lattice data, instead 
in the unconfined phase one finds the pattern 



21ogL 



(11) 



with a and b weakly dependent on T . This is shown in 
Fig. [1] for gluodynamics data [I3| 3 . The same pattern 
applies for the unquenched data of (3^|. For later refer- 
ence we quote the estimated continuum limit of the fitted 
values of a and b for gluodynamics [22| 



-0.19(5), 6=1.63(16). 



(12) 



3 An alternative approach is taken in |3l| to measure the renor- 
malized Polyakov loop in lattice QCD. 



A roughly constant a can be accommodated with per- 
turbation theory, and this is so quantitatively in the 
high temperature region. The second piece involves a 
power-like (as opposed to logarithmic) temperature de- 
pendence of non perturbative origin, since T c is propor- 
tional to Aqcd- On purely dimensional grounds it sug- 
gests a temperature independent non perturbative di- 
mension two condensate to be added in (A 2 a ). At zero 
temperature non perturbative condensates appear nat- 
urally within the operator product expansion approach 
[2lll33| through modifications in the propagators of parti- 
cles with respect to their perturbative counterparts. E.g . 
the standard dimension four gluon condensate [34], [35[ 
implies a term of the type (G^ t/ )/fc 6 to be added to the 
perturbative gluon propagator 1/fc 2 . Such non pertur- 
bative addition is irrelevant in the very high momentum 
regime but introduces sizeable power-like corrections as 
one approaches the strong interacting and non perturba- 
tive low momentum regime. A dimension two condensate 
has be en p roposed in the same zero temperature context 
in [36l . |37| |. This implies a gluon propagator of the type 
1/fe 2 + rrig/fc which, on dimensional grounds, yields a 
gluon exchange potential of the form Coulomb plus string 
tension, V ~ c/r + or (with a positive provided ttlq is 
also positive, corresponding to a tachyonic gluon squared 
mass -mg). As we will see immediately, such itiq gives 
rise to a dimension two condensate in the gluon field. 

Guided by this insight the same approach is taken in 
[22| to describe Polyakov loop data at not so high temper- 
atures above the deconfining transition. That is, a non 
perturbative piece driven by a dimension two operator is 
added in the finite temperature propagator: 



where 



D 00 (k) = D? (k) 



^o N o P (fc) 



m 2 D ) 2 



(13) 



(14) 



This expression nicely combines the effects coming from 
high frequency modes (the Debye mass from integration 
of the heavy thermal modes) and from the low frequency 
modes incorporated in the condensates. Thus we expect 
it to work in a window range above the phase transition. 
Using ([8]) at x = y, the new term produces 



NP 



(15) 



where the non perturbative condensate is related to ttVq 
through 



(Ka) 



NP 



l)m 2 G 



(16) 



8iririD 

This contribution is UV finite. When inserted in (0 it 
produces the term b(T c /T) 2 in the pattern Ijlljl . assum- 
ing that, ttiq or (A 2 a ) NP are temperature independent 
modulo radiative corrections. For future reference we in- 
troduce the notation 



C 2 



9 2 (Al a ) 



NP 



(17) 
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to denote the non perturbative dimension two conden- 
sate. Remarkably, the numerical value of the non per- 
turbative finite temperature condensate C2 turns out to 
be consistent [22| with independent determinations of the 
zero temperature condensate [38|, |39j, |40j . This is in favor 
of the consistency of the present approach. 



III. THE SINGLET QQ FREE ENERGY 

In this section we apply the previous model to the 
study of the free energy of a heavy quark-antiquark pair. 
The latter is obtained through the correlation of two 
Polyakov loops, 



-F eff (r,T)/T _ 



— trCl(x) — trftW 



(18) 



where r = \x — y\ is the separation between quark and an- 
tiquark. The free energy contains the finite temperature 
qq potential as well as the quarks selfenergies: 

F qq {r,T) = V qq (r,T) + 2E,(T), (19) 

where (twice) the selfenergy 

2E,(T) = lim F qq (r,T) (20) 

r— >oo 

is finite in the deconfincd phase and in this case the po- 
tential vanishes at large separations. In the confining 
phase the selfenergy would be infinite and F qq (r,T) is 
directly identified with the qq potential V qq (r,T). Of 
course, at large separations the two Polyakov loops in 
(I18p become uncorrelated implying that the quark (or 
antiquark) selfenergy is related to the expectation value 
of the Polyakov loop through 



L(T) 



-£„(T)/T 



(21) 



Because the interaction between the color singlets tr ft 
and trf^ requires the exchange of at least two gluons 
(each carrying a Debye screening factor eT mDT at large 
separations), this so called color averaged free energy is 
less suited to lattice computation than the singlet free 
energy 



-F 1 (r,T)/T 



i-tr (fl(x)tf(y)) 



(22) 



This quantity defined in the Coulomb gauge acquires a 
gauge invariant meaning (4~0 | and lattice calculations for 
it are available in the literature [42]. The correspond- 
ing qq potential is dominated by a single gluon exchange 
falloff at large separations. 

Following precisely the same approach described in the 
previous section for the Polyakov loop, we obtain the 
relation analogous to Q 



'^tr(n(x)nt(»)) 



(23) 



modulo 0(g 6 ) corrections. Once again, based on the 
cumulant expansion (30j these corrections vanish in the 
large N c limit. From the relation Fi = V\ + 2E 9 similar 
to (|19[) 4 we identify the two contributions 



Vi(r,T) = 



2iV c T 
2 



(A 0ia (x)A ,a(y)), 



AN r T 



(24) 



and in particular the results of the previous section for 
the Polyakov loop expectation value are recovered with 
(at LO) 



g 2 (N? 1) 

167r7V r 



-rh D T 



1 



2m n T 



(25) 



For the color singlet potential we derive the expression (at 
LO and neglecting any scale dependence in the coupling 
constant) 



Vi(r,T) 



9 2 {Nl 1 ) (1 
8ttN c 



1 



2rhrj T 



(26) 



In these results the Polyakov loop is already renormal- 
ized by dimensional regularization. In the corresponding 
lattice computation the bare Polyakov loop has to be 
renormalized using a suitable prescription. Specifically, 
in [l(J S3] the following criterion is used: the multiplica- 
tive renormalization of the Polyakov loop translates into 
an additive renormalization in the qq free energy, 



rcnorm 



(r > T) = F^™(r,T)+c(n 



(27) 



therefore, because at short distances the renormalized 
free energy ought to be temperature independent, a suit- 
able r-independent (but T-dependent) constant c(T) is 
added in such a way that at short distances the (renor- 
malized) finite temperature free energy F\ (r, T) has a 
maximum overlap with the zero temperature one F\ (r, 0) . 
The contact between the universal zero temperature 
curve and the temperature dependent one occurs for all 
separations below a distance of the order of the Debye 
length [10]. This device allows a quite accurate lattice 
determination of the selfenergy, and hence of the renor- 
malized Polyakov loop. In practical terms the above cri- 
terion means that, if the free energy is expanded in pow- 
ers of the separation, 



Fi(r, T) = v /r + Vl + v 2 r + 0(r 2 ) 



(28) 



(neglecting the logarithmic dependence from the running 
of the coupling constant), the constant term v\ should be 
absent (i.e. the Cornell prescription). This is because the 



Due to SU(A r c ) invariance, (Q) is a multiple of the identity matrix 
and so E 9 is the same quantity in F q q and F\. 
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dominant Coulombian term vq /r is temperature indepen- 
dent, and V\ vanishes in F\(r, 0) which is of the Coulomb 
plus string tension form. (A non-vanishing v\ would re- 
quire an inexistent "dimension one" condensate.) 

As it turns out, our formulas P5|) and pi?)) automat- 
ically satisfy this Cornell-Bielefeld prescription for the 
free energy. I.e., in F\ = V\ + 2T, q the term v\ is al- 
ready zero without any further renormalization. This is 
most remarkable since the selfenergy used a dimension- 
ally regularized definition of the UV divergent integral 
/ d 3 kDQ () (k) leading to (fTU)) . whereas the potential is 
UV finite for any non vanishing value of r. Actually, 
this was to be expected because, as is well-known, di- 
mensional regularization does not introduce any power- 
like dependence in addition to those already existing in 
the theory, so a dimensionful v± could not be introduced 
in the calculation. By the same token we expect that 
the Cornell-Bielefeld prescription will be automatically 
satisfied also at higher orders, i.e., beyond the Gaussian- 
like approximation, provided dimensional regularization 
is used. Note also that the fact that the Cornell-Bielefeld 
prescription used in lattice is satisfied by our approach 
justifies the comparison of our calculation of the renor- 
malized Polyakov loop, based on single loops, with that 
in lattice, based on the correlation of pairs of loops. 



IV. COMPARISON WITH LATTICE DATA 

A. Duality between the zero temperature 
QQ-potential and the Polyakov loop 

In this subsection we will analyze within our model 
the qq singlet free energy with respect to the two asymp- 
totic behaviors of zero temperature (at finite separation) 
and large separation (at finite temperature). Neglecting 
momentarily any scale dependence in both the coupling 
constant and the non perturbative condensate, we obtain 
(at LO) 



string tension 6 



Fx(oo,T) 



1 



-g m D T 



C 2 1 



8irN c * 2N C T 
for the singlet free energy at infinite separation, and 



(29) 



Fi(r,0) 



(AT C 2 - 1) g 2 1 



C 



2N r 4tt r 2N, 



Z + T^F^or (30) 



for the free energy at zero temperature (i.e., V qq (r)) 5 . 
The first thing to observe is that the latter takes the stan- 
dard zero temperature form of a perturbative Coulom- 
bian term plus a linearly raising confining term or, with 



5 Although the model is devised for the deconfined phase, nothing 
prevents us from taking the formal zero temperature limit in our 
expression for the free energy. 



2N, 



■C21TID 



(31) 



On the other hand, another obvious property in (1291) 
and (|30[) is the formal resemblance between the two 
asymptotic free energies with the identification r ~ 1/T. 
More precisely 



F^oo.T) =F 1 (r = l/m D ,0). 



(32) 



This relation indicates a remarkable and striking duality, 
within our model, between the quark selfenergy and the 
zero temperature qq potential. Such symmetry can be 
traced to the form of the propagator _Doo(^; T), namely, 
at least at LO its dependence on momentum and tem- 
perature comes in the combination k 2 + rn? D , regardless 
of the value of the coupling constant and the non pertur- 
bative condensate. Unfortunately the data are not avail- 
able in momentum space. This makes difficult to check 
the scaling of the two point function directly. However, 
when the momentum k is traded by the separation r, this 
symmetry can still be seen as a duality between the two 
asymptotic quantities Fi(po,T) and -Fi(r, 0). 

Of course the abovementioned duality between the 
roles played by k and mo will be disturbed by any scale 
dependence in the coupling constant and in the conden- 
sate corresponding to scaling violations. Such scale de- 
pendence implies a possible r and T dependence in those 
quantities and will be studied further below. Here we 
take g and C2 as constant (scale independent) but ad- 
mit for them the possibility of different values in the two 
asymptotic regimes (r,T = 0) and (r = 00, T). For the 
coupling constant this assumption is sustained by the 
data: the bosonic string model fits extremely well the 
zero temperature qq potential, with an error at the level 
of the 1% in the interval 0.15 fm < r < 0.8 fm [Hi ■ In 
this model the coupling constant takes a fixed value 



9(r,0) = g r 



7T 

2 



(33) 



Of course asymptotic freedom requires g(r, 0) to go to 
zero as r decreases, however this is only visible for sep- 
arations below 0.1 fm [Hj]. 7 On the other hand, in the 
large separation regime (r = 00, T), the Polyakov loop 
data favor a constant g (i.e. a constant a in (llljl ) rather 
than the perturbative running one. The constancy of 
a is unveiled only after the power correction has been 
subtracted off. Once again, asymptotic freedom requires 



At finite temperature, the 0(r) piece of the free energy is propor- 
tional to vOq — m? D indicating a weakening of the string tension. 
Also it confirms that the size of the contact between the zero and 
finite temperature curves is of the order of l/nip. 
Actually, the coupling g(fi) = tt/2 matching the string bosonic 
model corresponds to a NLO evolution scale 1/fi ~ 0.1 fm for 



A 



QCD 



240 MeV. 
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g(oo,T) to decrease as T increases but this perturba- 
tive regime occurs at very high temperatures only. This 
"freezing" of 5(00, T) to a constant value at intermediate 
temperatures is more thoroughly analyzed from the qq 
potential lattice data in the next subsection. So we take 



.9(00, T) = g T 
with gx obtained from the relation 

N? -la 2 . 



For N c = 3, N 



f = and a = —0.19 this gives 
g T = 1.21(11). 



(34) 



(35) 



(36) 



The fact that both g r and gx become frozen at intermedi- 
ates scales in the two regimes is in support of the duality 
discussed above. 

Let us now turn to the condensate, to see whether the 
relation (|3ip between the string tension and the conden- 
sate is verified. From the gluodynamics value of b in (fT2|) 
one obtains a condensate 



C 2 = (0.84(4) GeV) 2 



(37) 



When this is combined in (|3ip with g = g r (needed in 
171 d ) , it yields 



a = (0.43(2) GeV) 2 



(38) 



Remarkably, this is consistent with the accepted value of 
the string tension a as (0.42 GeV) 2 extracted, e.g., from 
the slope of the Regge trajectories 43] . In view of this, we 
will assume that the non perturbative condensate takes a 
common value in the two asymptotic regimes (r, T = 0) 
and (r = 00, T). 

Taking into account that g, but not the condensate, 
takes two different values in the two regimes, the duality 
relation (|3"2"|) gets modified since one has to use g equal 
to gr in (f2"9"| and equal to g r in ([3D)) . 8 This yields the 
refined duality relation 



Fi(oo,T) = 7 - 3 Fi(r = 7/m B ,0) 



(39) 



where 



7 



1/2 



m D = g T (N c /3 + Nf/Sf^T . (40) 



To see how this relation works, we plot in Fig.[5]lattice 
data points available for the two quantities: on the one 
hand twice the quark selfenergy, iq(oo,T) from [l(|, and 
on the other the zero temperature qq potential F\ (r, 0) 
from [TT| . To compare both quantities we use the re- 
lation r <-> 7/771,0 and the rescaling factor 7 3 indicated 
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FIG. 2: Verification of scaling relations in (I39|l . Black and 
white squares: 7 s Fi(oo,T) from lattice data in [10]. Circles: 
lattice data for the zero temperature qq potential from [TTj ] . 
Both quantities are plotted as functions of the separation r 
using the relation r <-> 7/mo. 



in p9p . 7 is computed using the numerical values of g r 
and gT quoted in (|33p and (|36p . which correspond to a 
continuum extrapolation of those parameters. 9 The re- 
markable agreement (modulo simple rescaling) between 
the two seemingly unrelated quantities suggests that this 
duality might actually be a feature of QCD, uncovered by 
our model. (Of course, any duality should break down at 
temperatures below the transition, since nothing is ex- 
pected to happen to the zero temperature qq potential 
for separations above 1/T C .) The same data are shown 
in Fig. [21 using the dimensionless quantities F\{oa,T)/T 
and (scaled) rF x (r,0), as functions of {T c /T) 2 . The 
straight line patterns (|29p and ([30]) are clearly exhibited, 
indicating the dominant role played by a non perturba- 
tive dimension two operator in the data. 



B. Singlet free energy at finite r and T 

In order to further check our model against non per- 
turbative information we will use the high quality lattice 
data of [43] , where the singlet qq free energy is computed 
for quenched QCD (iV c = 3) at various temperatures in 
the confining phase for temporal sizes N T = 4, 8. 

In order to make a comparison we use the following 
parameterization of the potential 



Vi(r,T) 



g 2 1 c 2 

3tt r 2N C T 



(41) 



Including the g in m^j, but of course not that in CV 



9 Note that the lattice data for Fi(r,0), unlike those of Fi(oo,T), 
have been extrapolated to the continuum. 
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FIG. 3: Same lattice data as in Fig. [2] for the dimensionless 
quantities Fi(pa,T)/T (squares) and BrV q q(r) (circles), as 
functions of (T c /T) 2 . (B — rriDj 4 .) Both the duality and 
the straight line pattern are exhibited by the data. 



FIG. 5: Points: values of <?t extracted from a fit to the 
Polyakov loop data [l^, 5H (circles) and from Fi (r, T) lat- 
tice data [43] (squares). LO formulas have been used and 
C 2 = (0.90(5) GeV) 2 . The curve displays the perturbative 
running coupling constant at three-loops. 
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FIG. 4: Points: lattice data, for iV T = 8, for the singlet free 
energy as a function of rT, at T/T c = 1.05, 1.5, 3, 6.06, 9, and 
12 (from top to bottom, respectively) from Ref. [43|. Curves: 
Fi(r,T) = Vi(r,T)+F 1 (oo,T) of our model. g(r) (to be used 
in Vi(r,T)) is extracted from the zero temperature potential 
[HI]. In addition C 2 = (0.9 GeV) 2 , and g T = 1.26. 



where rhp is related to g as in ([4]) with N c = 3, Nf — 
0. Furthermore, we take an r-dependent g consistent 
with the zero temperature gg-potential in namely, 
for each r, g(r) is taken so that 



9 2 (r) 
3nr 



(42) 



Following [llj, for r below 0.1 fm the perturbative value 
of V q q(r) to three-loops is used, and above 0.1 fm, g(r) — 
7r/2 from the bosonic string model. As noted before 
this hybrid model reproduces remarkably well the lat- 



tice data. The results are presented in Fig. 0] where the 
model is compared with N T — 8 lattice data (the results 
for iV r = 4 data are qualitatively similar). For clarity 
we represent Fi instead of Vy, and use rT as the inde- 
pendent variable. For the temperature dependent shift 
Fi(oo,T) the expression (|29|) is used with — 1.26, 
which corresponds to a value a = —0.21; this is the LO 
perturbative value around T — 6 T c . For the condensate 
we take C2 = (0.9 GeV) 2 which gives a fair account of 
Fi(oo,T) at lower temperatures. (Note that the points 
in the figure are the actual lattice data, without any con- 
tinuum extrapolation.) For consistency the same value 
of C2 is used in V\ (r, T) but the precise value is not very 
important in the potential (also it is not very important 
for F\(po, T) for high temperatures). As we can see from 
Fig. 21 the overall pattern of the lattice data is reasonably 
well reproduced by our model. Given the small error bars 
in the data, a more detailed agreement (small x 2 /pdof) 
was not to be expected without further sophistication of 
the model. As was noted in Section HH Fi(oo,T) (or 
equivalently, the Polyakov loop) cannot be described at 
all without introducing a dimension two condensate. The 
same holds here for the potential: the overall shape of 
lattice data for V\ (r, T) cannot be reproduced even ap- 
proximately if C2 is set to zero, moreover, the same value 
of the condensate extracted form F\{oo, T) works for the 
potential too. This feature is expected to survive in more 
elaborated models. 

We have considered also a fit to the data allowing a T 
dependent <?t (but keeping C2 constant). This is depicted 
in Fig. [5l The fit using all data for the free energy is of 
course more stable than that using only the Polyakov 
loop, although both sets yield consistent results. The er- 
ror bars are very small at high temperatures and come 
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from the lattice data and from allowing a variation range 
C 2 = (0.90±0.05 GeV) 2 in the condensate. We find an al- 
most flat dependence of gr with the temperature, much 
flatter than expected from a perturbative running. As 
noted, a nearly constant gx is in favor of the duality re- 
lation displayed in Fig. [2] and El since g r is also frozen 
in the bosonic string model, gr is somewhat larger than 
the perturbative value. It is noteworthy that a larger 
renormalization of the coupling constant, of about 50%, 
has also been observed in [42|, |45| if the observed Debye 
mass, extracted from the exponential falloff of the poten- 
tial, is compared to its perturbative value. In our model 
the exponential fall at large separations in controlled by 
g r rather than gr- Since g r = ir/2 is rather large, as 
compared to the perturbative value, this is qualitative 
agreement with the findings in [12, 5H (see also [IH), 
however, g r is somewhat smaller and does not run. 

The present analysis of the singlet qq free energy is car- 
ried out at LO, and further studies will have to take into 
account higher orders in the loop expansion. We certainly 
do not expect these higher order corrections to reproduce 
the effect of condensates. In addition, corrections com- 
ing from the condensate anomalous dimension are also 
expected to be relevant in more accurate descriptions; 
at the present stage this effect on the non-perturbative 
side should play a comparable role to loops corrections. 
Finally, higher dimensional condensates could also be 
present and affect the free energy. Unfortunately, as our 
previous analysis reveals [22| , the extraction of a dimen- 
sion four condensate from the lattice data of the Polyakov 
loop is not reliable at present. It would contribute to (TTTj) 
with a term (T c /T) 4 . 

V. CONCLUSIONS 

In the present paper we have investigated the role of 
dimension two condensates in the interrelation between 
the Polyakov loop at finite temperature above the phase 
transition and the static quark-antiquark potential at 
zero temperature. This extends our previous model for 
the renormalizcd Polyakov loop in the deconfincd phase 
which proved very successful and enabled a quantita- 
tive determination of the BRST invariant dimension two 



gluon condensate [22|. Here we have studied the conse- 
quences of the model when applied to the heavy quark- 
antiquark free energy in the unconfined phase of QCD. 

In our approach the perturbative and non-perturbative 
contributions are explicitly separated and an analytic and 
simple expression for the free energy which interpolates 
between the short and long distances regime is derived. 
The separation turns out to accommodate quite naturally 
a smooth running coupling when the model is confronted 
to lattice data, which speaks in favor of the consistency 
of the approach. 

We observe a natural consequence of the presence of 
the dimension two condensate if one restricts to LO in 
the perturbative contributions and effects associated to 
anomalous dimensions are disregarded. Although in gen- 
eral the free energy of a quark-antiquark pair depends 
both on the distance r and the temperature T, there 
are some regions where this double dependence presents 
strong correlations. Amazingly, the regimes of finite tem- 
perature and long distances (from which the Polyakov 
loop is defined) and the zero temperature and finite dis- 
tances (from which the zero temperature qq potential 
is obtained) have a striking dual functional dependence 
which is verified on a quantitative level on the lattice if 
the simple change of variable r «-> 1 /mj is assumed. The 
very fact that the non-perturbative contributions are en- 
coded in the dimension two condensate suggests that our 
result holds up to mild logarithmic radiative corrections 
and it would be interesting to quantify the observed small 
discrepancies in a more systematic manner. Despite all 
these caveats we believe that the quantitative success of 
the duality relations provides a further evidence on exis- 
tence and relevance of dimension two gluon condensates. 
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